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Executive  summary 


This  document  is  a  deliverable  for  the  study  Launching  of  Electromagnetic  Sur¬ 
face  Waves  on  Axial  Cylindrical  Reactive  Surfaces  with  Negative  Permittivity  under  the 
grant/cooperative  award  no.  FA8655-04-1-3005  and  was  issued  by  the  EOARD. 

In  this,  the  final  report,  the  excitation  of  the  dominant  TM  surface  wave  on  an  axial 
reactive  surface  composed  of  a  material  of  negative  permittivity  is  again  discussed  using 
a  different  and  more  accurate  analyses.  Here  we  tackle  the  problem  using  a  Wiener- 
Hopf  treatment.  The  geometry  is  slightly  modified  from  that  described  in  the  interim 
report  to  allow  variation  in  some  of  the  launcher  parameters.  The  surface  wave  launcher 
now  consists  of  a  metallic  cylindrical  surface  of  radius  b  coaxial  with  a  cylindrical  rod  of 
negative  permittivity  material  of  radius  a,  where  b  >  a.  The  reactive  surface  extends  from 
l  —  —oo  <  z  <  l,  and  the  conducting  surface  extends  from  — oo  <  z  <  0.  The  incident 
field  is  the  dominant  TM  mode  excited  in  the  coaxial  structure  and  propagating  in  the 
positive  2-direction. 

Numerical  results  are  presented  for  the  coupling  of  power  to  the  surface  wave,  the 
reflected  wave  and  the  radiated  field  on  the  plasma  rod.  The  calculations  were  performed 
over  a  wide  range  of  uip/u>  values.  The  case  of  waveguide  excitation,  i.e  b  =  a,  was 
compared  were  applicable  with  the  results  obtained  using  the  two  approximate  aperture 
distributions.  The  case  of  coaxial  excitation,  i.e.  b  >  a,  was  also  briefly  considered  using 
the  Wiener-Hopf  analyses.  The  calculated  launcher  efficiency  was  shown  to  increase  as 
the  b/a  ratio  increased. 
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1  Introduction 


1  Introduction 

1.1  Background 

The  problem  of  exciting  surface  waves  on  conventional  dielectric  rod  antennas  was  widely 
treated  in  the  1960  and  70’s  literature.  This  extensive  treatment  was  motivated  largely  by  the 
inherently  low  attenuation  and  large  bandwidth  available  in  some  surface  wave  launchers.  An 
excellent  survey  on  the  types  of  surface  waves  and  surface  wave  radiation  useful  in  dielectric 
rod  communication  antennas  has  been  collated  by  Chatterjee  in  [1,  2]  and  electrically  small 
dielectric  antenna  by  Luk  [3].  Recently,  researchers  have  returned  to  the  subject  with  a 
view  to  constructing  these  antenna  from  novel  dielectric  and  permeability  materials  that 
can  support  surface  waves.  These  materials  should  be  lightweight  and  robust  and  behave  as 
anisotropic  ferrite  and/or  plasma  analogues.  A  knowledge  of  the  properties  of  the  surface 
wave  and  plasma  material  properties  [4]  associated  with  these  structures  can  provide  useful 
insight  into  their  radiation  mechanism  and  characteristics. 

A  common  problem,  whether  using  a  conventional  dielectric  or  otherwise,  to  host  the 
surface  and  interior  fields,  is  the  efficient  excitation  or  launching  of  the  guided  electromagnetic 
surface  wave.  A  general  requirement  for  a  good  surface  wave  launcher  is  a  high  launcher 
efficiency  over  a  large  bandwidth.  As  was  pointed  out  in  the  interim  report,  when  the 
surface  wave  fields  are  of  infinite  extent,  the  launcher  in  theory  should  also  be  of  infinite 
extent  to  achieve  a  100%  launcher  efficiency.  Brown  demonstrated  theoretically  [5]  that  the 
launching  efficiency  of  a  finite-sized  launcher  could  approach  100%;  but  this  large  efficiency 
was  realised  only  at  the  expense  of  frequency  bandwidth. 

The  majority  of  reported  numerical  techniques  consider  the  excitation  source  to  be  either 
infinitesimally  short  electric  and/or  magnetic  dipole  current  sources.  A  finite-sized  aperture 
could  be  handled,  at  least  in  theory,  by  a  superposition  of  infinitesimal  sources.  In  practice, 
the  outcome  has  an  unknown  accuracy  since  the  formulation  fails  to  take  account  of  the 
losses  in  the  reflected  field  at  the  launcher.  Secondly,  because  the  aperture  distribution  may 
be  unknown.  In  the  interim  report  [6],  the  aperture  distribution  was  approximated  by  a 
“chopped”  surface  wave  distribution;  that  is,  the  fields  in  the  aperture  plane  are  assumed  to 
have  the  same  form  as  the  surface  wave  fields  within  the  aperture  and  are  assumed  to  vanish 
everywhere  outside  of  the  aperture.  Another  often  used  approximation  method  is  Kirchhoff ’s 
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approximation.  In  this  method  the  aperture  field  is  assumed  to  be  of  the  same  form  as  the 
unperturbed  incident  field.  This  distribution  was  also  computed  in  the  interim  report  [6] 
and  the  launcher  efficiency  results  were  found  to  be  poorer  than  those  obtained  using  the 
equivalent  aperture  technique.  The  principal  reason  for  the  discrepancy  was  attributed  to 
the  differing  radiation  patterns  obtained  as  the  ratio  uip/oj  was  increased.  The  radiated 
fields  were  computed  by  both  methods  by  an  integration  over  the  aperture  plane  (the  surface 
wave  modes  and  the  radiation  field  are  orthogonal  [7]).  The  accuracy  of  the  results  obtained 
using  these  approximation  techniques  is  usually  unknown  since  no  criterion  exists  which  can 
determine  the  extent  of  the  approximation. 

1.2  Scope 

The  purpose  of  this  report  is  to  (i)  present  an  exact  analysis,  using  Laplace  transform  and 
Weiner-Hopf  factorisation  techniques  (ii)  and  numerical  results  for  the  radiation  pattern, 
bandwidth  and  launching  characteristics  for  a  finite-sized  launcher.  These  results  will  then 
be  compared  with  (iii)  the  results  obtained  in  the  interim  report  [6]  for  the  “equivalent 
end-aperture”  and  Kirchhoff’s  approximations  to  determine  under  what  conditions  the  ap¬ 
proximate  techniques  are  valid. 


A  good  deal  of  the  ground  work  identifying  the  plasma  cylinder  dimensions,  useful  ma¬ 
terial  properties  and  resonant  eigenmodes  was  performed  in  [6].  The  analysis  considers  the 
geometry  shown  in  Fig.  1.  The  waveguide  launching  section  has  been  modified  from  that 
in  the  interim  report  to  allow  for  the  possibility  that  b  >  a  i.e.  coaxial  excitation  and  when 
b  =  a,  a  cylindrical  waveguide  excitation  as  in  the  interim  report.  This  modification  has  the 
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merit  of  generalising  the  analysis  and  design  parameters  but  comes  with  an  increased  cost 
in  analytic  complexity.  The  surface  wave  field  here  is  now  taken  account  of  by  specifying 
a  surface  impedance.  The  surface  impedance  is  the  ratio  of  tangential  electric  to  magnetic 
field  on  the  surface  of  the  wave  guiding  structure.  This  numerical  value  will  depend  on  the 
composition  of  the  structure  and  the  polarisation  of  the  propagating  field. 

The  structure  consists  of  a  reactive  negative  permittivity  cylindrical  surface  of  radius 
a  and  of  infinite  extent  in  the  ^-direction.  Coaxial  with  this  cylinder  is  a  thin,  perfectly 
conducting  surface  of  radius  b  for  z  <  0.  The  surface  wave  field,  radiation  field,  and  the 
reflected  field  will  be  computed  when  the  incident  field  is  the  dominant  TMoo  mode  in  the 
region  a  <  r  <  b,  z  <  0  propagating  to  the  right.  In  this  analysis  the  fields  in  the  region 
r  <  a  are  not  needed  and  if  they  exist  will  not  be  considered  in  the  calculations  of  power. 

2  Modal  Fields  and  Eigenfunctions 

2.1  Coaxial  Launcher  for  Axial  Cylindrical  Antenna 

This  problem  can  be  handled  with  the  Wiener-Hopf  technique.  The  procedure  to  be  followed 
is  similar  to  that  used  by  [8,  7].  The  analysis  will  be  carried  out  for  a  range  in  the  parameters 
of  the  structure  where  only  the  dominant  TMoo  mode  propagates  in  the  region  a  <  r  <  b,  z  < 
0.  The  total  field  in  this  case  can  be  derived  from  a  scalar  function  ip(r,z)  because  of  the 
circular  symmetry  of  the  structure.  The  function  ip(r,  z)  is  the  9  component  of  the  magnetic 
field.  In  the  region  a  <  r  <  b  and  z  <  0,  ip(r,  z)  can  be  expanded  in  a  series  of  eigenfunctions 
[9,  10], 


il>(r,z)  =  A0[Jx(-jpor)H^\-jpob) 

-M-jPob)Ji(-jp0r)H^2\-jp0r)\exp(-j0z) 

+  T%=iMJi(.Pnr)Hf\pnb) 

~MPnb)H[2\pnr)]exp(jnz)  (1) 

where  the  An  are  the  propagation  constants1.  The  Jo,Ji  represent  Bessel  functions  of 
the  first  and  second  kind  of  zero  and  first  order,  and  ,  Hq2'  are  Hankel  functions  of  first 
‘A  time  dependence  of  exp(juf)  is  assumed. 
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and  second  kind  or  order  zero.  The  presence  of  a  complex  argument  indicates  the  wave  is 
attenuated,  or  augmented,  in  the  outward  traveling  direction.  The  eigenvalues  pn  and  the 
propagation  constant  satisfy  the  equations 


2  2,72 

7o  =  Po  +  k0 

(2) 

1 n=pl~  fc0>for  n  >  0 

(3) 

_M-jPob)H^\-jp0a)  -  Jo(jp0a)H[2\-jp0b) 

Jpo ®  fn\  /o\  OLd 

J\(-jp0a)H $  \jp0b)  -  J0{jpob)H{  \jp0a) 

(4) 

Yo{Pnb)Y0(pna)  -  Jo(Pnb)Y0(pna) 

Poa  -r  /  x  Tf  '  .  —  aa,  for  n  >  0 

Jo{pnb)Yi(p0a)  -  Ji(pna)Y0(pna) 

(5) 

a  =  kXs/Zo 

(6) 

where  /cq  is  the  free-space  wave  number  and  Zq  the  characteristic  impedance  of  free-space 
and  To,  Yj  are  Bessel  functions  of  the  second  kind  or  order  zero  and  one.  The  eigenvalues  pn 
are  determined  by  the  boundary  conditions.  The  ratio  of  the  axial  component  of  the  electric 
field  to  the  azimuth  component  of  the  magnetic  field  [11]  must  equal  the  surface  reactance 
jXs  along  the  reactive  surface  at  r  =  a.  Along  the  surface  of  the  perfect  conductor  r  =  b  the 
axial  component  of  electric  field  must  vanish.  These  conditions  have  been  built  into  Eqn. 
1.  If  surface  reactance  tends  to  zero,  then  the  first  term  in  Eqn  1  will  reduce  to  the  TMoo 
TEM  like  wave  associated  with  a  coaxial  line  and  the  remaining  terms  will  reduce  to  the 
TM0n  modes  for  a  coaxial  line.  Alternatively,  if  b  =  a,  i.e.  cylindrical  waveguide  excitation 
then  the  first  term  in  Eqn.  1  reduces  to  zero  and  we  are  left  with  only  the  second  term 
again  permitting  only  the  T M0n  modes,  of  which  only  the  dominant  TMqi  is  considered.  It 
should  be  recalled  that  the  total  field  consists  entirely  of  TMqu  therefore  it  is  only  necessary 
to  specify  the  surface  reactance  for  one  polarisation  of  the  electric  field. 

2.2  Axial  Cylindrical  Antenna 

The  total  field  in  the  open  antenna  section  consists  of  discrete  modes  guided  by  the  reactive 
structure  plus  the  radiated  field.  Since  the  field  in  the  coaxial  portion  of  the  structure  is 
made  up  entirely  from  circularly  symmetric  TM  modes,  then  the  field  excited  in  the  open 
structure  must  also  be  circularly  symmetric,  the  only  mode  that  possess  circular  symmetry 
is  the  TMqi  mode. 
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The  total  field  in  the  open  structure  consists  of  two  components;  the  guided  TMq\  surface 
wave  plus  the  radiated  field.  The  surface  wave  must  be  of  the  form  [12], 


B0h[2)  (~jh0r)  =  exp(-j(30z)  (7) 

where  Bq  is  the  complex  amplitude  constant.  The  eigenvalues  ho  and  the  propagation 
constant  /3q  are  related  as, 


0o  — 


(8) 


The  eigenvalue  ho  satisfies  the  boundary  eigenfunction  equation  obtained  from  Eqn  6. 


jh0a 


H^(-jh0a) 
H[2\-jh0a) 


=  aa 


(9) 


The  resulting  eigenvalues  and  surface  reactance  values  are  consistent  with  those  obtain 
from  the  approximate  methods.  The  formal  solution  is  now  carried  out  using  the  Laplace 
transform  and  Wiener-Hopf  techniques  [11],  First,  the  total  field  ip(r,z)  is  decomposed  into 
an  incident  ipi(r,z)  and  scattered  field  ips(r,z)  components, 


ip(r,z)  =ipi(r,z)+ips(r,z)  (10) 

The  incident  field  from  the  coaxial  section  of  the  structure  will  be  of  the  form  TMq, 

z )  =  1^-{Ji{~jPor)H<0‘\-jpob)  -  M-jpob)H[2\-jp0r)}  exp(-j70z)  (11) 

The  scattered  field  beyond  z  >  0  satisfies  Helmholt’s  equation  expressed  in  cylindrical 
coordinates  by  Eqn.  12,  the  boundary  conditions  on  the  negative  permittivity  surface  and 
the  perfectly  conducting  surface  are  enforced  by  Eqns.  13  and  14,  respectively.  In  addition 
the  scattered  field  must  also  conform  to  the  radiation  condition,  which  requires  solutions  for 
the  scattered  field  to  be  divergent  waves  at  infinity. 
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S2(j>s  1  Scj)s 

Sr 2  r  Sr 


Sr2  +  (k0  1/r‘2)0s  =  0 

(12) 

~T~(r</>s)  +  acj)s  |r=:(1=  0 
r  or 

(13) 

~~  ~F~  (^^s)  |r=M<0=  0 
r  or 

(14) 

2.3  Laplace  Transformation  of  Scattered  Field 

The  scattered  field  can  be  expressed  as  a  bi-lateral  Laplace  transform  along  with  the  associ¬ 
ated  boundary  conditions, 


<t>(r,/3)  =  <f)+ (r ,  (3)  +  <j>  ( r,P )  (15) 

where, 

roo 

<j>+(r,p)  =  (f)(r,  z)  exp(— /?,  z)dz  (16) 

Jo 

<j)-(r,f3)=[  <p(r)z)exp(-p,z)dz  (17) 

J  —  OO 

The  scattered  field  can  then  be  recovered  by  the  inversion  integral, 

V>( r , z)  =  j  <f>{r,  P)  exp(/3,  z)dP  (18) 

To  make  (j>+{r,0 )  and  (p~  (r,  P)  analytic  functions  of  0  in  the  complex  p  region,  the  free 
space  wave  number  fco  will  be  made  complex.  This  is  equivalent  to  introducing  losses  into 
the  medium  surrounding  the  structure.  Let  jfco  =  jk'0  +  Tcq  ,  where  k'0  and  fco  are  real.  In  the 
final  solution  k'0  can  be  made  equal  to  zero  to  recover  the  results  for  the  lossless  case.  Since 
fco  is  complex,  the  propagation  constant  70  and  0q  will  be  complex.  Therefore,  in  general 
3lo  =  ji  +  31o  and  30 0  —  jPo  +  0o-  If  should  be  remembered  that  0'o ,  0$,  i'0  and  Jq  are  all 
real. 

Taking  the  Laplace  transform  of  Eqns.  12,  13  and  14  yields, 


Page  6  of  24 


Morrow 


2.4  Solution  of  Transformed  Scattered  Field 


+  “  +  (P2  +  k0  -  l/r2)<l>  =  0  (19) 

-  y~ ( r<t> )  +  acj)  |r=0=  0  (20) 

r  or 

r^(r^_)  +a^  lr=b=  °  (21) 

where  /3  is  restricted  to  the  range  —jg  <  Re(f3)  <  Jq. 

2.4  Solution  of  Transformed  Scattered  Field 

In  the  region  a  <  r  <  b  the  proper  solution  of  Eqn.  19  is  of  the  form, 


w  p)  =  Mh-  [A^2)(Aa  +  aH^{\a)]Jx{\r)  -  [AJ0(Aa)  +  aMXa)\H[2) (\r) 
[XH^2)(Xa  +  aJ?j2)(Aa)]  Ji(A6)  -  [AJ0(Aa)  +  <*Ji(Aa)]tf{2)(A6) 

The  solution  is  an  even  function  of  A  consequently  either  branch  of  A  can  be  selected.  In 
the  region  r  >  b  the  proper  solution  of  Eqn.  19  is, 

V>(r,/?)  =  ^(&+  (23) 

where  A  =  (fcg  +  /32)1/2.  We  will  now  determine  expressions  for  the  unknown  coefficients 
<j>{b~ ,P)  and  <j>(b+,/3)  in  Eqns.  22  and  23  by  using  the  discontinuity  condition  on  the  scattered 
field  at  r  =  b.  This  is  achieved  using  some  auxiliary  functions  J+(b,j3),  J~(b,f3)  and  E+(b,f3), 
where 


POO 

J+(b,P)=  /  [^s{b+,z) -ib(b~,z)]exp(-/3z)dz 

Jo 

(24) 

J~(b,P)=  [  [ips{b+,z) -ip(b~,z)]exp(-pz)dz 

J — OO 

(25) 

£+(M)  =  -f(n/>)  |r=6 

r  or 

(26) 

Is  should  be  mentioned  that  J~(b,fi)  is  the  Laplace  transform  of  the  electric  current  on 
the  outer  conductor  of  the  coaxial  guide  at  r  =  b  associated  with  the  scattered  field.  The 
function  E+(b,/3)  is  the  transform  of  the  axial  component  of  the  total  electric  field  evaluated 
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at  r  =  b.  The  function  J+(b,p)  can  be  deduced  for  the  region  Re(J3)  >  —  Jq  by  consideration 
of  the  boundary  conditions. 


J+(b,p) 


1 _ 

P  +  31o 


From  the  definitions  of  J+(b,/3)  and  J  (&,/?)  it  is  apparent  that 


(27) 


P(b+,P)  -  <t>(b-,P)  =  J+(b,f3)  -  J~(b,P) 
The  left  side  of  Eqn.  29  may  be  reduced  to, 


(28) 


<Kb+,p)-cl>(b~,p)=F(p)E+(b>p) 

where  F(0)  after  analytic  manipulation  is  explicitly  written  out  to  give, 


m  = 

nbX2  ffW(Xb) 


Hg(Aa 


(29) 


(30) 


X[l42\\a)J0{Xb)  -  H^\xb)J0{Xa))  +  a[H{2) {Xa)J0(Xb)  -  H^2)(Xb) Ji(Xa)] 

The  function  F(p)  can  be  shown  to  be  analytic  in  the  strip  —  7q  <  Re{P)  <  -Jq. 
Substituting  Eqn.  29  into  Eqn.  31  and  using  the  result  for  J+(b,f3)  from  Eqn.  28  gives, 


F(0)E+(b,p)  =  J~(b,p)  +  — —  (31) 

P  +  Jlo 

2.4.1  Wiener-Hopf  Factorisation 

Eqn.  31  can  be  solved  for  the  functions  E+{b,p)  and  J~(b,  /3)  by  performing  a  Wiener- 
Hopf  factorisation  on  the  function  F(/3);  i.e.  F{(3)  will  be  expressed  as  the  ratio  of  the  two 
functions  F+(P)  and  F~(0),  where  F+  is  analytic  and  nonzero  for  Re{fi)  >  jg  and  F~  is 
analytic  and  non-zero  for  Re{/3)  <  Jq.  Some  of  the  important  properties  of  the  functions  are 
listed  below; 

•  The  function  F+(/3 )  has  a  single  zero  at  0  =  —jp  and  a  branch  point  at  (3  =  ~jko- 
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•  The  function  F  (/ 3 )  has  an  infinite  number  of  zeros  in  the  complex  /3-plane  and  a  branch 
point  at  the  P  =  jko.  The  zeros  of  F~(/3)  are  located  at  (3  —  jyo  and  f3  =  7„(n  >  0). 

•  The  function  F+(/3);  is  of  the  order  (3~ at  infinity,  and  the  function  F~{(3 )  is  the 
order  /3-1/2. 


The  decomposition  of  the  function  F((3)  into  the  ratio  of  F+(/3)  to  F  {(3)  means  we  can 
transform  Eqn.  31  to, 


F+(p)E+(b,p)  =  } 

P  +  J'YO 


F-Q3)-F-(— j70) 
(3  +  j7o 


F-((3)J~(b,l3 ) 


(32) 


lm(P) 


|  /h  RVr 

J 

i 

^  I+fh 

.  ..  w, 

-yo” 

/ _ 

W 

yo” 

N 

F(P) 

/ 

Figure  2:  Complex  plane  diagram  illustrating  the  areas  were  the  various  transforms  are 
analytic. 

Fig.  2  shows  the  complex  /3-plane  and  illustrates  the  boundaries  of  the  various  analytic 
functions.  For  the  solution  of  the  scattered  field  to  be  unique  conditions  on  the  edge  or 
boundary  must  be  enforced  at  r  =  b  and  2  =  0.  These  conditions  require  the  axial  component 
of  the  electric  field  to  be  of  the  order  of  21/2  at  the  edge,  which  makes  the  transform  of  the 
electric  field  (i.e. ,  E+(b,(3))  of  the  order  of  /3-1/2  as  P  — >  00.  Similar  conditions  exist  for  the 
asymptotic  form  of  the  current  at  the  edge.  This  condition  requires  J~(b,P)  to  be  of  the 
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order  of  P  1  as  /3  — >  — oo.  Applying  Lioville’s  theorem  to  enforce  edge  conditions,  Eqn.  32 
can  be  equated  to  zero,  which  yields; 


£+(&,/?) 


F  {~ho) 

p+mp+ho) 


(33) 


2.4.2  Scattered  Field  Inversion  Integral 

The  coefficient  ip(b+,/3)  can  now  be  expressed  in  terms  of  the  known  function  E+(b,p)  by 
using  Eqns.  22,23,  and  26, 


W,/3) 


g|2)((fcg  +  /32)1/2b)F-(-j7 q) 

(fc02  +  P2)1/2H^({kl  +  0*)V*b)F+(f3)(P  +  no) 


(34) 


F  (~j 7o) 

F-(/3)(/3+no) 

The  scattered  magnetic  field,  tp(r ,  2),  may  now  be  computed  by  using  the  inversion  inte¬ 
gral, 


ip{r,z)  = -J—  [  <p{r,f3)exp((3z)dz  (35) 

Aft 3  Jc 

The  inversion  contour  C  is  shown  in  Fig.  3  and  must  be  located  in  the  strip  —/yg  < 
Re(/3)  <  7g  and  be  on  the  sheet  of  the  Riemann  surface  that  corresponds  to  the  choice 
Im(fcg  +/32)1/2  <  0  in  the  strip.  The  branch  cuts  are  selected  as  radial  straight  lines  starting 
from  the  branch  points  /3  =  ±jkQ.  For  the  lossless  case  the  branch  cuts  are  located  on  the 
imaginary  /3  axis. 

To  compute  the  discrete  spectrum,  the  amplitude  of  the  modes  must  be  evaluated  at  only 
one  radius.  It  seems  sensible  to  evaluate  the  field  at  the  radius  r  =  b  since  the  radial  varying 
factor  of  the  general  solutions  are  normalised  to  unity  at  r  =  b. 

Taking  the  region  a  <r  <b  and  z  <  0,  the  amplitude  of  the  scattered  field  at  r  —  b  is, 

ip(b,z)  = J^(b~,P)expci{Pz)dz  (36) 

where  ij;(b~,P)  is  given  by  equation  Eqn.  35.  The  contour  C  is  closed  in  the  right  half 
of  the  /3-plane  with  a  semi-circle  of  infinite  radius  deformed  around  the  branch  cut  as  shown 
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lm(P) 


Figure  3:  The  Inversion  Contour  in  the  Complex  Plane. 

in  Fig.  3.  The  details  of  the  integration  are  omitted  for  brevity,  but  are  based  upon  the 
technique  outlined  in  Chapter  4  of  [13],  for  integration  of  isolated  saddle  points.  Suffice  to 
say  the  integration  along  the  semi-circle  and  along  the  branch  cut  can  be  shown  to  be  zero 
for  z  <  0.  Thus,  the  integral  over  the  original  contour  must  equal  2irj  times  the  sum  of  the 
residues  of  the  enclosed  poles  in  the  right  half  /3  plane.  All  the  poles  of  the  ip(b~,/3 )  in  the 
right  half  are  due  to  the  zeros  of  the  F~((3).  Therefore, 


ips{b,z) 


F  (-j7o) 


2j7o 


sf-(0) 

m 


\p=. 


■  exp(j70z)  +  Y, 


F  (-j 7o) 


■3  7o 


(7 r 


1/3= 


■  exp(7„z)  (37) 


The  first  term  of  Eqn.  37  represents  the  dominant  TMqq  coaxial  mode  in  the  reflected 
field  and  the  summation  represents  the  evanescent  TMon  modes  in  the  reflected  field.  For 
z  <<  0  the  reflected  field  consists  entirely  of  the  first  term,  the  T Moo  mode. 

The  field  in  the  region  z  >  0  evaluated  at  r  =  b  can  be  found  by  using  the  same  inversion 
integral.  This  time  the  contour  is  closed  in  the  left  half  of  the  /?  plane  with  a  semi-circle 
of  infinite  radius  that  is  deformed  around  the  the  branch  cut.  The  integration  along  the 
semi-circle  can  be  shown  to  be  zero  for  z  >  0.  Therefore,  the  integral  must  equal  2nj  times 
the  sum  of  the  residues  of  the  enclosed  poles  in  the  left  half  /3  plane  plus  the  branch  cut 
integral.  The  poles  of  ip(b~,P)  in  the  left  half  0  plane  are  due  to  the  zeros  of  F+([3)  at 
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r  =  j(3 o  and  the  pole  appearing  explicitly  in  Eqn.  35  at  ft  =  -jjo-  Thus  for  2  >  0, 


ips{b,z)  = 


F  (~ho) 


-  jH[2\-jk0b ) 


(tf  -  PZy/*HP(-jhb)('y0  -  Po)^  1/3=100 


exp  (jpoz) 


(38) 


+  exp(— 7702)  +  branch  cut  integral 


The  first  term  in  Eqn.  39  represents  the  TM$\  surface  wave  mode  in  the  open  section 
of  the  antenna.  The  second  term  is  the  scattered  field  which  nulls  the  incident  field  in  the 
region  2  >  0. 

The  branch  cut  integral  or  continuous  part  of  the  spectrum  gives  rise  to  the  waves  as¬ 
sociated  with  the  radiated  field.  The  pattern  of  the  far  zone  radiated  field  is  evaluated 
approximately  by  the  saddle  point  method  of  integration,  the  contour  C  is  deformed  into  a 
contour  of  steepest  descent.  Transforming  r  and  2  from  cylindrical  to  a  spherical  co-ordinate 
frame  with  the  orientation  given  in  Fig.  1  the  expression  for  the  far-field  radiation  is, 


ips(p,4>)  = 


_ F  (— J7o)  exp  —j(kpp  +  7r/2) _ 

nF+(-jk  sin.  (p)(HQ2\kb cos  <j>))(jo  ~  k  sin  cj>)(kp  cos  <j>) 


(39) 


The  total  field  can  now  be  found  by  adding  the  scattered  field  to  the  incident  field  as 
given  in  Eqn.  10. 


3  Power  Flow  and  Antenna  Efficiency 

The  energy  transported  by  the  field  is  computed  using  the  Poynting  theorem.  The  only 
components  of  the  total  field  that  transport  energy  are  the  TMqo  incident  field  (Eqn.  11), 
the  TMqo  reflected  field  (first  term  of  37),  the  TMqi  surface  wave  field  (first  term  of  Eqn. 
39)  and  the  radiated  far  field,  Eqn.  39.  4>(r,z)  is  the  9  component  of  the  total  magnetic 
field,  the  electromagnetic  fields  f?,  7?  are  given  by, 


t  = 


jwe0 


-V  x  jjjae 


1}  =  Ipdg 


(40) 
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Thus  the  power  associated  with  the  total  field  is, 

~P  =  l/2Re^E!  x  if  -ds  (41) 

=  - - Im  [  (V  x  ipag)  x  ip*a$  ■  ds 

2uje0  Js 

where  ip*  denotes  the  complex  conjugate  and  s  is  a  surface  normal  to  the  surface  wave. 
For  the  case  of  the  incident,  reflected  and  surface  wave  fields,  the  only  component  of  the 
curl  of  ip  do  that  is  of  interest  is  the  radial  component  given  by  —Sip/5z.  The  power  guided, 
reflected  and  radiated  are  therefore, 


Pine  =  —  /  \ipi\2rdr 

(42) 

weo  Ja 

Pref  -  —  [  \ipr\2rdr 

(43) 

weo  Ja 

Psurf  =  —  [  \ips\2rdr 
weo  Ja 

(44) 

where  the  subscripts  refer  to  the  incident,  reflected  and  surface  wave  components,  respec¬ 
tively,  of  the  total  field.  For  the  case  of  the  radiated  field  the  component  of  ip  do  that  is  of 
interest  is  given  by  The  power  radiated  is  therefore, 


7T  ko 

Prad  =  -  /  \lprad\2p2  COS  (pd(p  (45) 

W£0  J-w/2 

As  defined  in  the  the  interim  report  the  quality  of  the  surface  wave  launcher  is  computed 
from  the  definition  of  launcher  efficiency  as  a  function  of  frequency.  The  launcher  efficiency 
was  defined  in  [6]  as  the  ratio  of  the  surface  wave  to  the  total  power  radiated  from  the 
aperture  of  the  launcher. 


V  laun 


Pflowi^  —  &) 


(46) 


Pflow{r  o)  T  Pflow{r  ■-> 

This  calculation  also  provides  an  opportunity  to  check  the  validity  of  the  two  approximate 
techniques  calculated  for  the  T Mqi  mode  against  the  results  of  an  exact  analysis. 
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4  Computed  Results 

Results  were  computed  for  the  case  of  the  cylindrical  waveguide  launcher  with  the  same 
dimensions  and  plasma  parameters  used  in  the  interim  report  (i.e.,  A  =  60.0cm.  /  =  3.0cm 
and  a  =  0.1cm).  To  replicate  the  excitation  of  the  plasma  rod  from  cylindrical  waveguide 
b/a  =  1. 

The  eigenvalues,  radiated  fields  and  power  transport  for  the  principal  TMoi  mode  were 
computed.  This  was  performed  on  the  Dell  Precision  360  Workstation  using  the  Matlab 
programming  language.  Matlab  possesses  several  very  useful  features;  the  Bessel  and  Hankel 
functions  used  were  “built  in”  and  conform  to  those  most  frequently  used  by  engineers  [14], 
it  automatically  selects  the  principal  value  of  an  analytic  function  and  lastly  has  a  toolbox 
that  can  handle  steepest  descent  integrations.  Table  1  shows  the  computed  permittivity  and 
surface  impedances  of  the  harmonic  cylindrical  fields  over  the  same  ratios  of  ljp/uj  used  in 
the  interim  report.  The  values  are  in  close  agreement  with  those  calculated  in  the  interim 
report. 


Plasma  Freq  Ratio 

(Wp/w) 

Permittivity 

(F/m) 

Surface  Impedance 

W 

5.4 

-28.1 

J71.1 

9.0 

-81.2 

J42.0 

100 

-10,000 

J3.7 

105 

-1.1  x  10n 

~  jo 

Table  1:  Table  of  Surface  Impedance  versus  Plasma  frequency  ratios  at  r  —  a  for  TMo\  mode 
on  the  plasma  rod  (A  =  60  cm,  l  =  3.0  cm,  a  =  0.1  cm). 


4.1  Radiated  Fields 

Figs.  4  and  5  show  the  radiated  fields  from  the  launcher  for  the  case  b/a  =  1  and  ojp/uj 
values  corresponding  to  those  given  in  Table  1.  (The  curves  have  been  normalised  to  the 
maximum  value  of  the  power  density  to  aid  comparison  with  those  in  the  interim  report). 
The  figure  illustrates  how  the  radiated  field  changes  for  different  values  of  surface  reactance  or 
equivalently  negative  permittivity.  The  field  pattern  shape  changes  with  increasing  negative 
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permittivity,  which  are  in  broad  agreement  with  those  predicted  in  the  interim  report  by 
the  equivalent  two-aperture  end  fire  analysis,  although  there  are  some  subtle  differences  that 
begin  to  emerge  when  uv/uj  >  102.  The  most  notable  difference  is  the  presence  of  a  small 
broadside  lobe  in  the  main  beam  when  ljp/uj  =  105.  It  appears  likely  that  another  mode, 
probably  the  HE\\ ,  is  being  partially  excited,  propagated  and  radiated;  and  some  of  this 
energy  has  been  picked  up  in  the  radiation  contour  integration.  The  inference  drawn  in  the 
interim  report  is  verified  here  that  the  approximate  equivalent  aperture  analysis  method  took 
more  precise  account  of  the  fields  present  than  the  Kirchhoff  method. 

4.2  Power  Propagated  and  Radiated 

The  results  for  the  powers  propagated  and  radiated  are  shown  in  Fig.  6,  for  a  value  of 
v/uj  =  0.01.  Where  v  is  the  collision  frequency  in  the  plasma  and  u>  the  radian  rf  excitation 
frequency2.  The  ratios  of  forward  traveling  (z  directed)  surface  wave  power,  reflected  power, 
and  radiated  power  to  the  incident  power  are  given  in  percentage  terms  as  a  function  of  up/u 
over  the  usual  range.  The  results  show  that  the  waveguide  launcher  is  efficient  in  launching 
the  TMoi  surface  wave  and  as  the  ratio  uip/u>  increases  the  reflected  power  decreases. 

4.3  Waveguide  Launcher  Efficiency 

Fig.  7  shows  the  launching  efficiency  vs  u>p/ui  for  the  structure  calculated  based  on  the 
“chopped”  surface  wave  distribution.  The  Wiener-Hopf  result  is  presented  along  with  the 
results  obtained  for  the  equivalent  end  aperture  and  Kirchhoff  surface  methods.  The  launch¬ 
ing  characteristics  both  exact  and  approximate  have  the  same  general  behavior  as  shown  in 
Fig.  7.  Closer  inspection  indicates  the  Wiener-Hopf  consistently  predicts  higher  launching 
efficiency  than  those  computed  using  the  approximate  methods.  This  feature  if  carried  over 
to  a  practical  antenna  design  would  deliver  less  losses.  With  regard  performance  bounds, 
the  approximate  methods  yield  quite  accurate  results  over  the  range  10_1  <  lup/uj  <  101. 
Beyond  this  point  the  Wiener-Hopf  results  diverge  from  the  approximate  and  this  is  thought 
to  be  commensurate  with  the  decrease  in  reflected  power  shown  in  Fig.  6. 

2  these  terms  and  their  relationship  to  the  plasma  frequency  up  and  permittivity  er  were  described  in  the 
interim  report. 
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End  Aperture  Radiation  Pattern 


(a) 


End  Aperture  Radiation  Pattern 


(b) 

Figure  4:  Radiation  pattern  due  to  fields  on  the  radial  surface  of  the  cylindrical  plasma 
supporting  the  TMqi  mode  for  various  u>p/ u  ratios.(a)  wp/w  =  5.4  (b)  lop/w  =  9.0. 
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End  Aperture  Radiation  Pattern 


(c) 

End  Aperture  Radiation  Pattern 


(d) 

Figure  5:  (Contn).  Radiation  pattern  due  to  fields  on  the  radial  surface  of  the  cylindrical 
plasma  supporting  the  TMqi  mode  for  various  up/w  ratios. (c)  u)p/oj  =  100  (d)  uip/ui  =  1  x  105. 
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Figure  6:  Ratio  of  Powers  Coupled  from  the  Incident  wave  to  the  forward  travelling  Surface 
Wave,  Reflected  Surface  Wave  and  Radiated  Field  for  the  TMq\  mode. 

4.4  Coaxial  Launcher  Efficiency 

In  studying  a  specific  design  of  launcher,  i.e.  the  waveguide  launcher,  we  were  able  to 
determine  that  a  modestly  efficient,  electrically  small  plasma  antenna  was  practically  feasible. 
However  other  forms  of  launcher  do  exist  which  might  deliver  higher  coupled  powers  to  the 
antenna.  Perhaps,  the  most  effective  of  which  is  a  coaxial  launcher.  Since  by  varying  the 
inner  to  outer  ratio  of  the  radius  of  the  plasma  rod  to  metallic  cylinder  we  can  vary  the 
impedance  matching  between  the  propagating  and  radiating  wave-fields  (as  suggested  in 
Section  1.2). 

The  outer  radius  b  was  increased  to  a  ratio  of  b/a  =  2.3.  The  ratios  of  2.3  corresponds  to 
a  coaxial  characteristic  impedances  of  50f2  and  is  not  intended  to  be  an  optimum  impedance 
match  over  the  whole  range  of  plasma  ratios.  It  is,  however,  intended  to  provide  some 
impedance  matching  over  the  critical  region  of  the  antennas  operation  i.e.,  101  <  uip/u>  < 
102.  Fig.  8  compares  the  results  for  the  launching  efficiency  of  the  waveguide  and  coaxial 
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Launcher  Efficiency 


Figure  7:  Computed  Launcher  Efficiency  for  Cylindrical  Waveguide  Excitation  with  the 
T Moi  mode  using  the  End  aperture,  Kirchhoff  and  the  Wiener-Hopf  methods. 

launchers.  From  the  figure  it  is  clear  their  is  no  significant  improvement  below  ujp/ui  <  1, 
beyond  this  region  the  launching  efficiency  is  observed  to  increase  slightly. 

5  Conclusions 

An  antenna  composed  of  a  negative  permittivity  material  and  excited  via  a  waveguide 
launcher  has  been  investigated  mathematically.  Two  independent  treatments  of  the  problem 
have  been  given;  In  the  interim  report  two  approximate  analyses  using  the  Kirchhoff  and 
equivalent  end  aperture  method  were  applied  and  in  this  report  an  exact  method,  using 
the  Wiener-Hopf  technique  developed.  The  exact  solution  provides  results  with  bounded 
accuracy  and  are  considered  more  rigorous. 

•  The  launcher  efficiency  was  computed  for  all  three  methods  based  upon  the  “chopped 
surface  wave”  definition.  On  the  approximate  methods,  the  equivalent  end  aperture 
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Launcher  Efficiency 


Figure  8:  Computed  Launcher  Efficiency  for  Cylindrical  Waveguide  and  Coaxial  Excitation 
(waveguide  b/a  =  1,  coaxial  b/a  =  2.3). 

results  were  more  accurate  than  the  Kirchhoff  solution  when  compared  with  the  Wiener- 
Hopf  solution. 

•  Both  Wiener-Hopf  and  end  aperture  analyses  indicate  the  possibility  for  constructing  an 
electrically  small  antenna  from  a  negative  permittivity  material  that  is  fairly  radiation 
efficient.  The  analyses  here  further  endorses  the  summary  conclusion  made  in  the 
interim  report. 

•  The  Wiener-Hopf  solutions  consistently  calculated  slightly  higher  launching  efficiencies 
up  to  (jjp/u  =  101  and  beyond  this  point  did  not  cut  off.  One  reason  suggested  for 
this  behaviour  is  when  ujp/ui  >  101,  the  internal  fields  are  being  gradually  occluded 
from  the  interior  plasma  and  energy  transport  is  occurring  only  on  the  surface,  were 
the  impedance  losses  are  low.  In  the  approximate  methods  losses  were  assumed  and 
computed  throughout  the  cross-section  of  the  interior  guide. 
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•  For  completeness  a  coaxial  launcher  design  was  also  investigated.  This  was  demon¬ 
strated  to  provide  more  effective  matching  of  the  launcher  to  the  plasma  geometry. 
The  results  indicate  that  this  could  enhance  launcher  efficiency  by  another  10%  for 
uip/uj  >  1.  As  the  plasma  ratio  increases  the  radiation  patterns  followed  the  waveguide 
behaviour. 


6  Recommendations  for  Future  Work 

The  following  recommendations  are  made; 

•  It  would  seem  appropriate  to  experimentally  substantiate  some  of  the  claims  made  in 
the  conclusions.  This  could  be  achieved  via  another  cooperative  research  program. 
Reported  experiments  with  negative  permittivity  antennas  have  concentrated  on  free- 
space  plasmas  [15,  16].  However,  these  antenna  suffers  from  high  losses  and  present 
many  practical  implementation  issues.  The  idea  proposed  here  is  to  construct  a  passive 
plasma  analogue  material3,  excite  it  using  the  coaxial  launcher  and  characterise  its 
performance.  The  intended  operating  band  would  be  in  the  20-1000  MHz  frequency 
spectrum  were  compact,  lightweight  wideband  aerials  for  mobile  communications  are 
often  required.  Transmit  power  would  be  restricted  to  a  few  10’s  of  Watts.  A  single 
antenna  is  suggested  were  the  engineering  novelty  should  be  to  provide  an  antenna 
much  more  compact  than  a  conventional  A/4  or  A/2  whip  aerial. 

•  It  is  unlikely  that  a  single  antenna  will  be  suitable  to  cover  the  whole  frequency  spec¬ 
trum  mentioned  above.  Therefore,  I  should  wish  to  extent  the  research  effort  to  either 
compact  clustered  or  active  plasma  analogue  variants.  This  latter  effort  would  address 
the  challenge  of  achieving  more  bandwidth  and/or  more  gain.  AFRL  should  advise 
if  such  a  study  is  relevant  to  their  current  requirements  and  how  best  to  initiate  the 
study. 

•  One  other  attendant  research  issue  is  the  antenna  operation.  The  radiation  from  any 
electrically  small  HF  antenna,  which  is  dipole  like,  will  be  effected  by  the  presence  of 

3  Some  interesting  papers  were  highlighted  in  the  Interim  Report  that  point  to  useful  passive  negative 
permittivity  materials  that  may  be  suitable  for  antenna  application. 
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the  equipment  operator,  the  ground  and  surrounding  environment.  It  would  also  be 
appropriate  to  investigate  experimentally,  or  otherwise,  the  effect  on  radiation  pattern 
and  antenna  system  losses  in  the  vicinity  of  dielectric  and  conductive  objects. 

•  To  make  appropriate  conference  and  journal  publications  based  on  this  work  and  citing 
the  EOARD  sponsorship.  It  is  noted  there  has  been  a  recent  call  for  a  special  issue  of 
the  IEE  on  metamaterials  and  novel  antennas. 
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